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ABSTRACT 

In cold dark matter cosmologies, small mass halos outnumber larger mass halos at any redshift. 
However, the lower bound for the mass of a galaxy is unknown, as are the typical luminosity of the 
smallest galaxies and their numbers in the universe. The answers depend on the extent to which 
star formation in the first population of small mass halos may be suppressed by radiative feedback 
loops operating over cosmological distance scales. If early populations of dwarf galaxies did form in 
significant number, their relics should be found today in the Local Group. These relics have been 
termed "fossils of the first galaxies" . This paper is a review that summarizes our ongoing efforts to 
simulate and identify these fossils around the Milky Way and Andromeda. 

It is widely believed that reionization of the intergalactic medium would have stopped star formation 
in the fossils of the first galaxies. Thus, they should be among the oldest objects in the Universe. 
However, here we dispute this idea and discuss a physical mechanism whereby relatively recent episodes 
of gas accretion and star formation would be produced in some fossils of the first galaxies. We argue 
that fossils may be characterized either by a single old population of stars or by a bimodal star 
formation history. We also propose that the same mechanism could turn small mass dark halos 
formed before reionization into gas-rich but starless "dark galaxies" . 

We believe that current observational data supports the thesis that a fraction of the new ultra- faint 
dwarfs recently discovered in the Local Group are in fact fossils of the first galaxies. 
Subject headings: cosmology: theory — galaxies: formation — stars: formation 



1. INTRODUCTION 

There are many questions that remain open in cosmol- 
ogy with regard to the mass, number and properties of 
the smallest galaxies in the universe. Have we already 
discovered the smallest galaxies in the universe or we are 
still missing an elusive but large population of ultra- faint 
dwarf galaxies? 

In cold dark matter (CDM) cosmologies most of the 
dark halos that formed before reion ization had masses 
small er than 10 8 - 10 9 M (e.g., iGnedin fc Ostrikerl 
1997). The small mass halos that survived tidal destruc- 
tion to the modern epoch, were they able to form stars, 
would constitute a sub-population of dwarf satellites or- 
biting larger halos. Small mass dark halos significantly 
outnumber more massive galaxies like the Milky Way and 
can be located in the voids between luminous galaxies 
fe.q- lHoeft eTaL|[200llRicottil [20091 ). 

However, until recently (i.e., before 2005) observations 
did not show a large number of satellites around massive 
galaxies like the Milky Way and Andromeda. This be- 
came known as the "missing galactic s atellite problem, 
(|Klypin et al.lll999t [Moore et all 11999ft . The voids be- 
twe en bright galaxies appear to be devoid of dwarf galax - 
ies (jKarachentsev et all 120041 12001 iTullv et alj [2006ft . 
While the abundance of dwarfs in large voids may 
not pose a probl e m to CDM cosmology, as shown by 
iTinker fc Conrovl (|2009ft . it is unclear whether the pre- 
dictions of the number of faint dwarfs in the Local Group 
is consistent with both the number of observed Milky 
Way dwarf satellites and the number of relatively iso- 
lated dwarfs in the local voids. 

Historically, the discrepancy between observation and 
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theory on the number of dwarf galaxies has been inter- 
preted in two ways: 1) as a problem with the CDM 
paradigm that could be solved by a modification of the 
dark matter properties - for instance b y introducing 
warm dark matter (e.o. JBode et al.ll2001ft - or 2) as an 
indicator of feedback processes that are exceptionally ef- 
ficient in preventing star format ion in small mass hal os, 
which remain mostly dark (e.g.. Haiman et al . 2000a). 

The recently discovered populat i on of ultra-faint 
dwarfs (iBelokurov et all I2006L 120071: llrwin et all 120071: 
Walsh et all 120071: iWillman et all l2005bUal: iZucker et all 
2006b lal: llbata et alJ l2007t iMajewski et al.l l2007t 
Martin et al.ll2006ft in combination with a proper treat- 



ment o f observational i ncomp le teness (iSimon fe Gehal 
20071 IKoposov et alJ [2001 iTollerud et alJ 120081 : 



Geha et al .1 12008ft has increased the estimated num- 
ber of Milky Way satellites to a level that can 
be more easily reconciled with theoretical expecta- 
tions. For instance, the suppression of dwarf galaxy 
formation due to intergalactic medium (IG M ) re- 
heating during reionization (iBabul fc Reed 119921: 



IShapiro et all 11994 
Thoul fc Weinberg! 1199a iQuinn et alJ 1 1996 
Navarro fc Steinmetzl 11997 
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Diikstra et all 120041: IShapiro et al] 12004 iHoeft et all 



200a lOkamoto et al 1120081 lRicottil l2009) , in conjunction 
with a strong suppression of star formation in small mass 
pre-reionization dwarfs, may be sufficient to explain the 
observed number of Milky Way satellites. In the near 
future we can hope to answer perhaps a more interesting 
question: what is the minimum mass that a galaxy can 
have? This is a non trivial and fundamental question in 
cosmology. Answering it requires a better understanding 
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of the feedback mechanisms that regulate the formation 
of the first galaxies before reionization and the details of 
the process of reionization feedback itself. 

The formation of the first dwarf galaxies - before 
reionization - is self-regulated on cosmological distance 
scales. This means that the fate of small mass halos (i.e., 
whether they remain dark or form stars) depends on local 
and global feedback effects. This type of galaxy feedback 
differs from the more familiar model operating in normal 
galaxies (e.g., SN explosions, AGN feedback, etc), where 
the feedback is responsible for regulating the star forma- 
tion rate within the galaxy itself but does not impact 
star formation in other distant galaxies. Rather, before 
reionization, each proto-dwarf galaxy reacts to the exis- 
tence of all the others. Different theoretical assumptions 
and models for the cosmological self-regulation mecha- 
nisms will, of course, produce different predictions for 
t he number and luminosity of the first dwarf galaxies 
(lHaiman et alj|2000al iRicotti et alJl2002U iWise fc Abell 
2008 UO'Shea fc Normanll2008t iRicotti et al.ll2008al) ~ 

We now introduce the basic concepts on how feedback- 
regulated galaxy formation operates in the early universe 
(i.e., before reionization): 

A cooling mechanism for the gas is required in order to 
initiate star formation in dark halos. In proto-galaxies 
that form after reionization this is initially provided by 
hydrogen Lyman-alpha emission. This cooling is efficient 
at gas temperatures of 20,000 K but becomes negligible 
below T ~ 10, 000 K. Later, as the temperature drops be- 
low 10,000 K, the cooling is typically provided by metal 
line cooling. In the first galaxies, however, both these 
cooling mechanisms may be absent. This is because the 
first dwarf galaxies differ when compared to present-day 
galaxies in two respects: 1) they lack important coolants 
- such as carbon and oxygen - because the gas is nearly 
primordial in composition, and 2) due to the smaller typ- 
ical masses of the first dark halos, the gas initially has 
a temperature that is too low to cool by Lyman-alpha 
emission. 

The gas in small mass halos with circular veloc- 
ity w V ir = (GAftotAvir) 1 / 2 ~ 20 km s _1 , where r vir is 
the virial radius - roughly corresponding to a mass 
Mtot ~ 10 8 M Q at the typical redshift of virialization - 
has a temperature at virialization T 10, 000 K. Hence, 
if the gas has primordial composition it is unable to 
cool by Lyman-alpha emission and initiate star forma- 
tion unless it can form a sufficient amount of primor- 
dial H2 (an abundance xh 2 ~ 10~ 4 is required). Because 
molecular hydrogen is easily destroyed by far-ultraviolet 
(FUV) radiation in the Lyman- Werner bands (11.3 < 
hv < 13.6 eV) emitted by the first stars, it is widely be- 
lieved that the maj ority of galaxies with y v - lr < 20 km s _1 
remain dark (e.g., [Haiman et aLll2000al ). However, sev- 
eral studies show that even if the FUV radiation back- 
ground is strong, a small amount of H2 can always form, 
particularly in relatively massive halos wi th virial tem- 
perature of several thousands of degrees (|Wise fc Abell 
l2007tlO'Shea fc Normanll200l . Thus, negative feedback 
from FUV radiation may only delay star formation in the 
most massi ve pre-reionization dwarfs r ather than fully 
suppress it ijMachacek et al.l l2001. 2003). We will argue 
later that hydrogen ionizing radiation (hv > 13.6 eV) in 
the extreme-ultraviolet (EUV) plays a far more impor- 



tant role in regulating the formation of the first galaxies 
than FUV radiation. Thus, in our opinion, models that 
do not include 3D radiative transfer of H and He ioniz- 
ing radiation cannot capture the most relevant feedback 
mechanis m that regulates gal axy formation in the early 
universe (|Ricotti et al.ll2002al |bl). 

After reionization, the formation of dwarf galaxies with 
tVir < 20 km s _1 is strongly inhibited by the increase 
in the Jeans mass in the IGM. Thus, according to this 
model, reionization feedback and negative feedback due 
to H2 photo-dissociation by the FUV background (im- 
portant before reionization) determine the mass of the 
smallest galactic building blocks. The resultant circu- 
lar velocity of the smallest galactic building blocks is 
iVi r ~ 20 km s _1 , roughly corresponding to masses 
Mtot ~ 10 s — 10 9 M . If this is what really happens in 
the early universe, the "missing Galactic satellite prob- 
lem" can be considered qualitatively solved because the 
predicted number of Milky Way satellites with w vir > 20 
km s _1 is already comparable to the estimated number 
of observed satellites after applying completeness correc- 
tions (although this model may still have problems re- 
producing the observations in detail). 

However, as briefly mentioned above, we have argued 
for some time that most simulations of the first galax- 
ies cannot capture the main feedback mechanism oper- 
ating in the early universe because they do not include 
a key physical ingredient: radiative transfer of H and He 
ionizing radiation. Our simulations of the formation of 
the first galaxies are to date the only simulations of a 
cosmologically representative volume of the universe (at 
z ~ 10) that inclu de 3D radiative trans fe r of H and He 
ionizing radiation (|Ricotti et ah 2002alfbl . I2008H) . Fig- 
ure [1] shows the evolution of ionized bubbles around the 
first galaxies in a cubic volume of 1.5 Mpc in size at 
rcdshifts z = 21.2, 17.2, 15.7, 13.3 from one of our simu- 
lations. The results suggest that negative feedback from 
the FUV background is not the dominant feedback mech- 
anism that regulates galaxy formation before reioniza- 
tion. Rather, "posit ive feedback" on H2 f o rmation from 
ionizin g radiation (| Haiman et al.l Il996at IRicotti et al.l 
l2001al ) dominates over the negative feedback of H 2 disso- 
ciating radiation. Hence, a strong suppression of galaxy 
formation in halos with v v - 1T < 20 km s^ 1 does not take 
place. In this latter case, some galactic satellites would 
be the fossil remnants of the first galaxies. Comparisons 
of simulated pre-reionization fossils to dwarf spheroidals 
in the Local Group show rema rkable agreement in prop- 
erties (IRicotti fc Gnedinll2005l hereafter RG05). Based 
on the results of the simulations, we also suggested the 
existence of the ultra-faint population before it was dis- 
cover ed about a year later (see RG05, iBovill fc Ricottll 
l200flh . 

1.1. Definition of "pre-reionization fossils" 

Throughout this paper we define "pre-reionization fos- 
sils" as the dwarfs hosted in halos with a maximum cir- 
cular velocity remaining below 20 km s _1 at all times 
during their evolution: v m&x (t) < 20 km s _1 . It will 
become clear in this paper that this definition is not di- 
rectly related to the ability of fossils to retain gas and 
form stars after reionization: in §2] we describe a mecha- 
nism in which small mass halos with v max (t) < 20 km s _1 
are able to have a late phase of gas accretion and possibly 
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Fig. 1. — 3D rendering of cosmological H II regions (fully ionized gas is red and partially ionized gas is blue) around the first galaxies 
in a bo x of 1.5 Mpc. The four boxes show a time sequence at redshifts z = 21.2, 17.2, 15.7, 13.3 for the simulation S2 from IHicotti et all 
(2002b). The rendering shows several tens of small size H II regions around the first galaxies (there are a few hundreds of galaxies in this 
volume). A movie of the same simulation shows that the H II regions are short lived: they form, expand to a size comparable to the large 
scale filamentary structure of the IGM and recombine, promoting the formation of molecular hydrogen inside the relic H II regions. 



star formation. 

Our definition of fossil reflects the special cooling mech- 
anisms and feedback processes that regulate star forma- 
tion and the number of luminous halos with t> m ax(£) < 
20 km s _1 , before and after reionization. In proto-fossil 
galaxies - even adopting the most conservative assump- 
tion of maximum efficiency of shock heating of the gas 
during virialization - the gas is heated to a tempera- 
ture below T ~ 10, 000 K. Thus the gas cannot cool by 
Lyman-alpha emission, a very efficient coolant. The cool- 
ing of the gas is dominated either by H2 roto- vibrational 
line emission or by metal cooling, im portant if the metal- 
licity exceeds Z ~ 10~ 3 Z w (e.g., [Bromm et al.l l200lt 
Santoro fc ShuTl[200l IRicotti et "al.M2008bt iSmith et al l 
2009ft. These coolants are much less efficient than 
Lyman-alpha emission. Moreover, H2 abundance and 
cooling is modulated and often suppressed by the FUV 
and EUV radiation fields. The FUV radiation in the 
H2 Lyman- Werner bands and hard ultraviolet radia- 
tion have large mean free paths with respect to the 
typical distances between galaxies, thus their feedback 
is global in nature. Qualitatively this explains why 



the first galaxies have low luminosities and low surface- 
brightness, similar to dwarf sp heroidal (dSph) galax- 
ies in the Local Group (RG05, iBovill fc Ricottil l2009t 
iSalvadori fc Ferraral [20090 . 

Simulations also show that stars in the the first galax- 
ies do not form in a disk but in a spheroid (RG05, 
IRicotti et al.ll2008bl ). A thin galactic disk is not formed 
because of the high merger rates and the low masses 
of dark halos in the early universe. Roughly, pre- 
reionization fossils have a mass at virialization M to t < 
10 8 M Q , assuming they form at z v \ r ~ 10, but their 
mass may increase by up to one order of magnitude by 
z = due to secondary infall (RGS02a, b. IRicotti et alJ 
l2008a| ). Secondary infall does not affect v maxi which re- 
mains roughly constant after virialization. 

1.2. Pre-reionization fossils and reionization 

The critical value of w max ,crit for which dwarf galaxy 
formation is suppressed by reionization feedback is close 
to the 20 km s -1 value that defines a fossil, but it is 
not necessarily the same value. Indeed, it can be sig- 
nificantly larger than 20 km s _1 if the IGM is heated 
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to T > 10,000 K (jRicotti et al.l [2000h . Thus, we ex- 
pect that the virialization of new "pre-reionization fos- 
sils" is strongly suppressed after reionization due to IGM 
reheating (i.e., they mostly form before reionization). 
However, pre-reionization fossils and dark halos with 
« mal < 20 km s _1 that virialized before reionization may 
accrete gas and, in certain cases, form new sta rs after 
reionization at redshifts z < 1 — 2 (|Ricottill2009| ). 

Unfortunately, the value of iWx.crit is uncertain due 
to our poor understandin g of the thermal history of the 
IGM (jRicotti et al.l 120001 ). The uncertainty surround- 
ing the IGM equation of state may partially explain the 
differences found in literature on the values for v max ,crit 
and the different levels of suppression of star formation 
as a function of the halo mass after reionization (e.g., 
[ Weinberg etU] Il997l : ICnedTnl l2lwlal IHoeft et al.ll200l 
IQkamoto et al.ll2008l ) . Regardless of assumptions for the 
reionization feedback model, one should bear in mind 
that no halo with i> max < 20 km s _1 can form stars af- 
ter reionization unless the gas in those halos has been 
significan tly pre-enriched with m etals. For instance, the 
model bv iKoposov et al . (2009) assumes star formation 
after reionization in halos as small as w max ~ 10 km s _1 . 
With this assumption they find that their model is con- 
sistent with observations of ultra-faint dwarfs but claim 
that fossils are not needed to explain the data. However, 
star formation in such small halos can only take place 
in a gas that was pre-enriched with metals, suggesting 
the existence of older populations of stars in those halos. 
Indeed, according to our definition, the smallest post- 
reioniz ation dwarfs with 10 k m s _1 < i> ma x < 20 km s _1 
in the IKoposov et al.l (|2009f ) model are "fossils". As 
stated above, fossils may also be able to form stars af- 
ter reionization due to a late phase of cold gas accretion 
from the IGM (|Ricottill2009[ ). 

1.3. Identification of pre-reionization fossils in 
observations 

Pre-reionization fossils are not easily identifiable be- 
cause i> max cannot be measured directly from obser- 
vations. Understanding the star formation history of 
dwarf galaxies may help in this respect, as fossils likely 
show some degree of suppression of their star forma- 
tion rate occurring about 12.5 Gyr ago due to reioniza- 
tion. However, their identification based on their star 
formation history may be complicated because some pre- 
reionization fossils in the last 10 Gyr may have had a 
late phase of gas accretion and star formation. The 
caveat is that star formation histories cannot be mea- 
sured with accuracy better than to within 1-2 Gyr and 
the accuracy becomes increasingly poorer for old stel- 
lar populations. Thus, it is impossible to prove whether 
an old population of stars formed before reionization 
(which happened about 1 Gyr after the Big Bang) or at 
z ~ 3, when the Milky Way was assembled. Nevertheless, 
ultra-faint dwarfs that show some degree of bimodality 
in their star formation history are candidates for being 
pre-reionization fossils. 

Acc ording to results by RG05 and iBovill &; Ricottil 
(2009), Willman 1, Bootes II, Segue 1 and Segue 2 do 
not lie on the luminosity-surface brightness relationship 
of simulated pre-reionization fossils. This result is based 
on the assumption that fossil properties are not modi- 



fied by tides. Their surface brightness is larger than the 
model predictions for objects with such low-luminosity. 
An as yet undiscovered population of ultra-faints with 
lower surface brightness is instead predicted by our sim- 
ulations. It is likely that the properties of the lowest 
luminosity ultra-faints may have been modified by tidal 
forces due to their proximity to the Milky Way disk. 

Although it is difficult to identify individual fossils, 
statistical arguments suggest that at least some ultra- 
faint dwarf galaxies are pre-reionization fossils. This is 
because the number of satellites from N-body simula- 
tions with v max (t) > 20 km s _1 is substantially smaller 
than the estimated number of observed satellites after 
completeness corrections. Admittedly the current theo- 
retical and observational uncertainties on the number of 
satellites are still large. However, if the estimated num- 
ber (after completeness corrections) of ultra-faint dwarfs 
increases further, the existence of pre-reionization fossils 
will be inescapably proven. This is especially the case if a 
population of ultra-faint dwarfs with luminosities similar 
to Willman 1, Bootes II, Segue 1 and Segue 2 but sur- 
face brightness below the current sensitivity limit of the 
SDSS - as predicted by our simulations - is discovered. 

The possibility of identifying the fossils of the first 
galaxies in our own backyard is very exciting. It would 
greatly improve our understanding of the physics in- 
volved in self-regulating the formation of the first galax- 
ies before reionization. Clearly, even the launch of the 
James Webb Space Telescope (JWST), would not yield 
the wealth of observational data on the formation of the 
first galaxies that could be obtained by studying ultra- 
faint galaxies in the Local Group. 

The rest of the paper is organized as follows. In §[2] we 
briefly review and discuss observational data on Galactic 
satellites, in §[3] we summarize the results of simulations 
of the formation of the first galaxies in a cosmological 
volume and the effect of reionization feedback on galaxy 
formation. In § [5] we discuss a recently proposed model 
for "late gas accretion" from the IGM onto small mass 
halos. In § [5] we compare the theoretical properties of 
simulated pre-reionization fossils to observations. In §[6] 
we compare different ideas for the origin of classical and 
ultra-faint dwarf spheroidals. We present our conclusions 
in§El 

2. OBSERVATIONS 

2.1. The ultra-faint satellites of the Milky Way and 
Andromeda 

Prior to 2005, the number of observed dwarf sa tellites 
of th e Milky Way and Andromeda was about 30 (|Mateol 
1998). One of the most evident properties of the dwarfs 
in the Local Group is a type segregation, with "gas free" 
dwarf spheroidal (dSph) galaxies distributed near the 
center of their host galaxy and gas rich dwarf Irregulars 
(dlrr) at larger distances from the galactic centers. No- 
table exceptions are the Magellanic Clouds that are dlrr 
less than 100 kpc from the center of the Milky Way and a 
few isolated dSphs like Tucana and Antlia. One popular 
explanation for this segregation is the transformation of 
dlrr into dSph due to tidal and ram press ure stripping 
as dwarfs fall toward the Milk y Way center (|Maver et al.l 
2007: Pcilarru bia et al.l 120081 ). In addition, simulations 
showed that the number of dark matter satellites of the 
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Milky Way with mass > 10 8 M© (i.e., with mass suffi- 
ciently large to expect star formation in them) was an 
order of magnit ude larger than the number of know n 
dwarf satellites (jKlvpin et all 119991 : iMoore et all 11999( 1 . 
This posed a problem for CDM cosmologies. 

Since 2005-2006 the number of known Local Group 
satellites has begun to increase dramatically, with 
the discovery of a new population of ultra-faint 
dwarfs. The new galaxies have been discovered by 
data mining the SDSS and other surveys of the halo 
around M31, resulting in the discovery of 14 new 
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1 and 11 new companions for M31 


(Ilbata et al. 




2007 


; IMaiewski et all 120071 IMartin et all 



120061 I2009H . Unofficial reports from members of the 
SDSS collaboration state that there are actually at least 
17 new ultra-faint Milky Way dwarfs, but several of them 
are as yet unpublished (anonymous referee's private com- 
munication) . 

The new Milky Way satellites have been slowly 
discovered since SDSS Data Release 2, with the 
most recent, Segue 2 discovered in Data release 7 



Adelman-McCarthv et al l 120061 12001. IKoposqv et"aL1 
2008H and iWalsh et alj (|2009( 1 systematically searched 
Data Releases 5 and 6, respectively. Due to the par- 
tial sky coverage of the SDSS, an d assuming isotropic 
distribution of sate llites (but see IKroupa et alj 120051 : 
IZentner et al.l I2005H . the total number of ultra- faint 
dwarfs in the Milky Way should be at least 5.15 time s 
larger than the observed number (| Toller ud et al.l [2008) . 
With this simple but conservative correction, the num- 
ber of Milky Way satellites within 400 kpc is about 
12 + 5.15 x (14 ± 3.7) ~ 84 ± 19. The quoted uncer- 
tainty is simply Poisson error due to the relatively small 
number of known ultra-faint dwarfs. 

In estimating the completeness correction for the num- 
ber of Milky Way dwarfs, one should account for se- 
lection effects inherent in the method used to find the 
ultra-faints in the SDSS data. In addition to complete- 
ness corrections for the survey's coverage of the sky, 
the most important selection effect is the total number 
of stars from the object seen in the survey, i.e., hori- 
zontal branch stars or main sequence and/or red giant 
stars for the lowest luminosity ultra-faints like Coma or 
Segue 1. This sets a limiting surface brig htness cutoff at 
rough ly 30 ma g arcsec" 2 for the SDSS (Koposo v et al.l 
I2008D (but see lMartin et all (|2008( 1 that find a limiting 
surface brightness about 6.4 times larger). There is also 
a distance-dependent absolute magnitude cutoff. The ef- 
ficiency of finding ultra-faint dwarfs by data mining the 
SDSS typically drops rapidly at Galactocentric distances 
beyond 50 — 150 k pc for the ultra- faints (depending on 
their luminosity) ([ Simon fc Gehal [20071: IKoposqv et all 



I27M iTollerud et al.l 12008b IWalsh et alll2009f) . Of the 
new Milky Way dwarfs, only Leo T is well beyond this 
distance threshold 1 , and 11 of the 14 new Milky Way 
satellites are within 200 kpc. 
The luminous satellites can be radially biased. So, the 

1 Leo T was found because it contains a young stellar population 
and gas. Otherwise, it would not have been identified as an ultra- 
faint dwarf due to its large Galactocentric distance. 



abundance of the faintest satellites within 50 kpc, that is 
the most complete sample, may not be used to correct for 
incompleteness at larger distances from the Galactic cen- 
ter without prior knowledge of this bias. And, of course, 
satellites of different luminosity and surface brightness 
will have different completeness limits. These selection 
biases have been considered in a paper by ITollerud et al.l 
( 2008) . Their study finds that there may be between 300 
to 600 luminous satellites within 400 kpc. Their estimate 
for the number of luminous satellites within a Galacto- 
centric distance of 200 k pc is between 176 to 330. 

Re c ent surveys of M31 ((Martin et alj|2006t Ilbata et all 
2007; Mar tin et al.ll2009f ) have covered approximately a 
quarter of the space around the M31 spiral. The sur- 
veys have found 11 new M31 satellites, bringing the to- 
tal number to 19. If we make a simple correction for the 
covered area of the survey, the estimated number of M31 
satellites, including the new dwarfs, increases from 8 to 
52 ± 13. 

The sensitivity limits of the surveys for Andromeda do 
not allow the detection of ultra-faint dwarfs that would 
correspond with those with the lowest luminosity found 
in the Milky Way. However, despite the fact that An- 
dromeda and the Milky Way are thought to have ap- 
proximately the same mass (within a factor of two), their 
satellite systems show interesting differences for even the 
satellites at the brighter end of the luminosity function. 
For instance, there are small differences in the galacto- 
centric distance distribution of satellites and in the mor- 
phology of the satellites (e.g., number of dlrr, dE and 
dSphs). 

2.2. Peculiar ultra-faint dwarfs 

Many of the newly discovered dwarfs are dSphs with 
a dominant old population of stars and virtually no gas, 
which makes them candidates for being pre-reionization 
fossils. However, there are notable exceptions that we 
discuss below that may not perfectly fit the properties 
of simulated "fossils". For instance, the dwarf galaxy 
Leo T resembles all the other ultra-faints but has gas and 
recen t star formation ((Irwin et al.ll2007l ISimon k Geha 
12007( 1. We argue that Leo T could be a true "fossil" 
with v max (t) < 20 km s , but may have experienced a 
late phase of gas accretion and star formation due to the 
mechanism discussed in § |H 

Leo T has a stellar velocity dispersion of (JLeoT — 
7.5 ± 1.6km s _1 jS imon fc Gehal [2001 . or an estimated 
dynamical mass of 10 7 M© within the stellar spheroid (al- 
though its total halo mass may be much larger). Leo T 
shows no sign of re cent tidal destructio n by either the 
Milky Way or M31 (|de Jong et all 12008( 1 and is located 
in the outskirts of the Milky Way at a Galactocentric 
distance of 400 kpc. Leo T photometric properties are 
consistent with those of pre-reionization fossils. On the 
other hand, the halo of Leo T could be sufficiently mas- 
sive (« max (t) > 20 km s _1 ) to retain or accrete mass after 
rcionization and not be a pre-reionization fossil. How- 
ever, as discussed in § 0] it is also possible that Leo T 
is a pre-reionization fossil that has been able to accrete 
gas from the IGM at late times due to an increase in the 
concentration of its dark ha lo and a decrease in the IGM 
temperature (|Ricot ti 2009). Under this scenario, Leo T 
stopped forming stars after reionization, but was able 
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to start accreting gas again from the IGM at z ^ 1 — 2 
and therefore would have a bimodal stellar population. 
Ide Jong et al.l (|2008f ) have found evidence for bimodal 
star formation in Leo T. Our model would explain why 
Leo T does not resemble other dlrr and is similar to 
dSphs and pre-reionization fossils, while not suffering sig- 
nificant tidal stripping. 

Willman 1, Bootes II, Seguel and Segue2 are among 
the lowest luminosity ultra-faint dwarfs discovered so far; 
however, they do not fit the typical properties of pre- 
reionization fossils (see § [5]). For instance Willman 1 
has a dynamical mass within the largest stellar orbit 
(r ~ 100 pc) of 5 x 10 5 M Q and a mas s-to-light ratio 
~ 470 , similar to other ultra- faint dwarfs (Willman et al. 
I2005a|) . However, given its low luminosity, Willman 1 has 
central surface brightness that is too large when com- 
pared to simulated "fossils". Simulated fossils with lu- 
minosities Ly < 10 3 L Q should have a typical surface 
brightness that falls below the detectio n limit of ~30 
mag arcsec" 2 estimated for the SDSS (|Koposov et al.l 
2008). Hence, the lowest luminosity fossils may still be 
undiscovered. Although the nature of the lowest lumi- 
nosity ultra-faints is unknown, due to selection effects 
they can only be found within ~ 50 kpc of the Galactic 
center. Thus it is possible, and perhaps to be expected, 
that their properties have been affected by tidal forces 
(|Penarrubia et alj 120081 : IMartin et al.ll2008f i. 

3. FORMATION OF FIRST GALAXIES IN CDM 

The first episodes of star formation in the universe 
are thought to take place at redshift z ~ 30 — 50, 
in the center of dark matter halos with typical mass 
Mtot ~ 10 5 - 10 6 M Q . The gas in these halos is 
metal free and simulations show that a single or bi- 
nary massive star per halo is formed ( Bromm et al.l 
200alAbel et al.ll2002l: ISaieo et al.ll2004l: iGaq et al.ll2005l: 
O'Shea fc Normanll2007HYoshida et al.ll2008t iTurk et all 
20091: iStacv et all 120091) . Such stars are called Pop III 
and their mass, although not well constrained by the 
simulations, is quite large: in the range between 20 M 
to a few 100 M Q . Whether or not we can refer to miniha- 
los containing a single massive star (or a binary star) in 
their center as the "first galaxies" is debatable. However, 
the crucial point to be made here is that there is a gap 
of 2 to 3 orders of magnitude between the typical halo 
mass in which Pop III stars are born (10 5 — 10 6 M Q ) and 
the typical mass of the population of dwarf galaxies that 
are not strongly influenced by radiative and reionization 
feedback (10 8 — 10 9 M Q ). The primordial dwarfs that 
fill the gap are those that we refer to as pre-reionization 
fossils. 

If the formation of pre-reionization fossils is not drasti- 
cally suppressed by radiative feedback, their number may 
be several orders of magnitude larger than the number of 
more massive dwarfs. This is because in CDM cosmolo- 
gies the number of dark halos per unit comoving volume 
roughly scales with the mass as N oc M7r. 

It is widely believed that nearly all halos with mass 
M to t > 10 8 — 10 9 Mq host luminous galaxies, although 
there can be substantial disagreement among theorists 
on their luminosity. However, most of the theoretical 
controversy rests in understanding the fate of the halos 
with mass between 10 6 — 10 8 M Q and the dominant feed- 
back that determines whether they become luminous or 



remain dark. We will elaborate on this statement in the 
next sections. 

3.1. Radiative feedback 

Simulating the formation of the first stars is a rela- 
tively well defined initial condition problem, given the 
cosmological parameters. However, these simple initial 
conditions must soon be modified to take into account 
the effects of other newborn stars, whose properties are 
still quite uncertain. The physics becomes more complex 
as competing feedback effects determine the fate of the 
first galaxies: radiative feedback regulates the formation 
and destruction of H2 and metals are injected into the 
IGM and into protogalaxies. 

3.1.1. Negative feedback from H2 photo-dissociating 
radiation 

The net effect of radiative feedback on the global star 
formation history of the universe before the redshift 
of reionization is uncertain. An FUV background (at 
energies between 11.34 eV and 13.6 eV) destroys H2, 
the primary coolant at the start of galaxy formation. 
The FUV radiation emitted by the first few Pop III 
stars is sufficient to suppress or delay galaxy forma- 
tion in halos with circular velocities v v - n - < 20 km s _1 
that are too small to cool by Lyman-alpha emission 
(IHaiman et alJ2000btlCiardi et alJ200ttlMachacek et all 
l2000tlO'Shea k Normanll2008D . Hence, according to this 
scenario, most halos with masses < 10 8 — 10 9 M Q re- 
main dark. More work is needed to quantify the level of 
suppression of galaxy formation and examine how these 
models compare to observations of Milky Way satellites. 

Figure [2] illustrates the effect of H2 dissociating radi- 
ation on the I GM. The two p anels s how slices through 
a simulation in iRicotti et~aL1 |2002b) at z = 19.44 and 
z = 18.5. The top- right tiles in the two panels show 
H2 abundance. At z = 19.4, the H2 has its relic abun- 
dance everywhere in the IGM except inside the dissoci- 
ation spheres around the first galaxies, where it is de- 
stroyed. At z = 18.5, the dissociation spheres are still 
visible, but the UV background starts to dissociate H2 
everywhere in the IGM except the denser filaments. 

3.1.2. Positive Feedback Regions 

Our main criticism for the "negative feedback" model 
is that it does not ta ke into account the eff e ct of hydrogen 
ionizing radiation (|Ricotti et al.l l2001bt iWhalen et all 
l2007h that, according to simulations, may indeed play 
a dominant role in regulating g alaxy formation before 
reionization (jRicotti et al. 2002a.b). Simulations includ- 
ing 3D radiation transfer show that star formation in the 
first small mass halos is inefficient, partially due to winds 
produced by internal UV sources. This produces galax- 
ies that are extremely faint and have very low surface 
brightnesses. However, our simulations show that a large 
number of ultra-faint dwarfs (a few hundred galaxies per 
comoving Mpc 3 ) form before reionization at z ~ 7 — 10. 
Hence, according to this model, the Local Group may 
contain thousands of ultra-faint dwarf galaxies. 

Ionizing radiation from the first stars enhances the 
production of H2 (we refer to this as "positive feed- 
back") by creating free electrons and promoting the 
formation of H _ , the main catalyst for the formation 
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Fig. 2. — The two (2x2) panels show slices through the most massive object in the simulation 64L05p2 in RGS02b at z = 19.4 and 18.5. 
The box size is Lf, ox = 0.7 comoving Mpc. Each one of the 2x2 panels shows in log-scale: the neutral hydrogen fraction (top left), the 
molecular fraction (top right), the gas overdensity (bottom left), and the gas temperature (bottom right). The sequence illustrates the 
evolution of a H2 dissociation sphere around a single source (panel at z = 19.4) and the dominance of the H2 dissociating background at 
2 = 18.5. 



of H 2 in a low metall i city gas (IShapiro fc Kand IT9871; 
Haiman Rees. fe Loebl Il996ri IRicotti. Gnedin. fc Shulll 
2001b; I Alvarez. Bromm. fc Shapiro! 120061 : iCiardi et alJ 
2006h . iRicotti et all (|2001b[ ) found that shells of H 2 



can be continuously created in precursors around the 
Stromgren spheres produced by ionizing sources and, for 
a bursting mode of star formation, inside recombining 
H II regions. We refer to these shells as "positive feed- 
back regions". This is because the catalyst H~, and 
hence H2 , is formed most efficiently in regions where the 
gas ionization fraction is about 50%. This local "posi- 
tive feedback" is difficult to incorporate into cosmological 
simulations because the implementation of spatially in- 
homogeneous, time-dependent radiative transfer is com- 
putationally expensive. 

Fig. [3] shows "positive feedback regions" in one of our 
simulations. The figure shows a slice through a simu- 
lation at 4 different times (at z — 17.3, 12.2, 11.3, and 
10.2). We recognize the two main processes that create 
H2 in the filaments. In the top- left frame at z = 17.3 we 
can see a "positive feedback region" as an irregular shell 
of H2 surrounding the H II region that is barely inter- 
sected by the slice. In the bottom-left frame (z = 11.3) 
two H II regions are clearly visible. Inside the H II re- 
gions, the H2 is destroyed. In the bottom-right frame 
(z = 10.2) the H II regions are recombining (demonstrat- 
ing that the star formation is bursting) and new H2 is 
being reformed inside the relic H II regions. A finer in- 
spection 2 of the time evolution of this slice shows that 
at least five H II regions form and recombine between 
z = 20 and z = 10 in this small region of the simulation. 

There are two reasons why our results are still contro- 
versial. First, our simulations do not yet have sufficient 
resolution to ensure their convergence. Second, there are 
no other published simulations to compare our results 
with. Only recently have some groups started to in- 
clude the effect of 3D radiative transfer on hydrodynam- 

2 Movies of 2D slices and 3D rendering of the sim- 
ulations are publicly available on the web at the URL: 
http: / /astro. umd.edu/~ricotti/movies. html 



ics (e.g., I Wise fc Abel 120071 12001 . However, currently 
there are no other simulations of the formation of the first 
galaxies in a cosmological volume suited for comparison 
with observati ons of dwarfs in the L ocal Group other 
than our own (IRicotti et all l2002allbl; IRicotti fc Gnedinl 
20051 iRicottfetal . 2008a; Bovill fc R icotti 2009). Hence, 
our results may differ from other numerical studies be- 
cause of the inclusion in the code of the effects of "pos- 
itive feedback regions" and galactic winds from ionizing 
radiation. 

Simulations by iWise fc Abel (|2008f ) include a self- 
consistent treatment of hydro and 3D radiative transfer 
that is more accurate than our approximate, but faster 
method. However, because the authors use ray-tracing 
for the radiative transfer, only a few sources of radia- 
tion can be simulated at the same time. This limits 
the volume and number of galaxies that can be simu- 
lated. Due to these limitations the simulations are not 
suited for comparison between the primordial dwarf pop- 
ulations and the ultra- faint dwarfs. In addition, at the 
moment, the aforementioned simulations do not include 
metal cooling and the formation of normal stars (other 
than Pop III). 

3.1.3. The simulations 

The simulation used for comparison to observations of 
ultra-faint dwarfs has been thoroughly described in Ri- 
cotti, Gnedin & Shull (2002a, 2002b) as run "256Llp3" . 
Here we remind the reader that the simulation includes 
256 3 dark matter particles, an equal number of baryonic 
cells, and more than 700,000 stellar particles in a box of 
size ~ 1.5 comoving Mpc. The mass of the dark matter 
particles in our simulation is 4930 M Q , and real comoving 
spatial resolution (twice the Plummer softening length) 
is 150 pc (which corresponds to a physical scale of 24 
parsecs at z = 8.3). This resolution allows us to resolve 
cores of all simulated galaxies that would correspond to 
the observed Local Group dwarfs. The stellar masses are 
always smaller than the initial baryon mass in each cell 
but can vary from ~ 0.6 h _1 M Q to 600 h _1 M with 



Fig. 3. — Same as in Fig. [2] except for a zoomed region of 0.125 2 h~ 2 Mpc 2 around the most massive object in the 64L05p3 simulation in 
RGS02b. In this time sequence of images (top: z = 17.3, 12.2 from left to right; bottom: z = 11.3, 10.2 from left to right) we recognize the 
two main processes that create H2 in the filaments: "positive feedback regions" in front of H II regions and the reformation of H2 inside 
relic H II regions. The bursting mode of the star formation is evident from the continuous formation and recombination of the H II regions 
in the time sequence of the slices. 



a mean of 6 h _1 M Q . Stellar particles do not represent 
individual stars but, in general, a collection of stars (e.g., 
OB associations). 

The simulation includes most of the relevant 
physics, including time-dependent spatially-variable ra- 
diative transfer usin g the OTVET approximation 
(jGnedin fc Abel 120011 ). detailed radiative transfer in 
Lyman-Werner bands, non-equilibrium ionization bal- 
ance, etc. In addition to primordial chemistry and 3D ra- 
diative transfer, the simulations include a sub-grid recipe 
for star formation, metal production by SNe and metal 
cooling. The code also includes mechanical feedback by 
SN explosions. However, we found that for a Salpeter 
IMF, the effect of SNe is not dominant when compared 
to fee dback produced by io nizing radiation from massive 
stars (jRicotti et all l2008al hereafter RGS08). The ef- 
fect of SN explosions is somewhat model dependent and 
uncertain because it is treated using a sub-grid recipe. 
Hence, the simulation analyzed in this work includes 
metal pollution but not mechanical feedback by SNe. 

In RG05, we included the effect of reionization in the 



simulation 256Llp3. Because the size of the simulation 
box has been fixed at ~ 1.5 comoving Mpc, the simula- 
tion volume is too small to model the process of cosmo- 
logical reionization with sufficient accuracy. We there- 
fore assume that the simulation volume is located inside 
an H II region of a bright galaxy at a higher redshift. 
Specifically, we introduce a source of ionizing radiation 
within the computational box, properly biased, which 
corresponds to a star-forming galaxy with the constant 
star formation rate of 1 solar mass per year (similar to 
star formation rates of observed Lyman Break Galaxies 
at z ~ 4, Steidel et al. 1999). The source is switched on 
at z = 9.0, and by z = 8.3 the whole simulation box is 
completely ionized. 

3.1.4. Summary of main results 

T he main results of th e simulations are the following 
(see lRicotti et al.ll2008al for details): 

1. Negative feedback. H2 photo-dissociation from 
FUV radiation, the main negative feedback 
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Fig. 4. — (Left). Fraction of baryons retained by each galaxy (normalized to the cosmic mean value, M bar = /j^ x M dm , »ucic Jbar 
^b/n m = 0.136) as function of the halo mass M^ m for run SI in RGS08 at z = 10 . The size of the dots is proportional to the fraction of 
stars f t = M*/M baI in each halo: from the largest to the smallest dots we have /* > 10%, 1 < /* < 10%, 0.1 < /* < 1% and /» < 0.1%, 
respectively. The plot illustrates the role of internal and external sources of ionizing radiation in reducing the gas retained by small mass 
halos. Luminous sources (with larger /„) retain less gas than dark halos due to winds driven by internal sources of radiation. (Right). 
Average star formation efficiency (/*) as a function of halo mass at z = 10.2 (left panel) and 12.5 (right panel) for run SI. We divide all 
halos into three groups: those at distance d < 8 kpc from the nearest luminous halo (solid histograms), those with 8 kpc < d < 50 kpc 
(dashed histogram), and those with d > 50 kpc (long-dashed histograms). The dot-dashed curve shows the fraction of luminous halos 
flum(Wd m ) as a function of the halo mass. 



thought to suppress the formation of the first galax- 
ies, is not the dominant feedback. If we modify the 
spectrum of the sources of radiation to artificially 
increase or decrease the FUV flux emitted by up to 
one order of magnitude, we do not find any appre- 
ciable effect on the global star formation history. 

2. Self -regulation. Feedback by hydrogen ionizing ra- 
diation (EUV) plays the key role. Galactic out- 
flows, produced by UV photo-heating from massive 
stars, and H2 formation/photo-dissociation induces 
a bursting star formation mode in the first galaxies 
that acts as the catalyst for H2 re-formation inside 
relic (recombining) H II regions and in the "precur- 
sors" of cosmological Stromgren spheres (i.e., posi- 
tive feedback regions). As a result, the formation of 
the first galaxies is self-regulated on a cosmological 
distance scale. It is significantly reduced by radia- 
tive feedback but it is not completely suppressed, 
even in halos with w max ~ 5 — 10 km s _1 . Note that 
our sub-grid recipe for star formation assumes a 
Schmidt law with a tunable efficiency parameter e* 
(the fraction of gas converted into stars per cross- 
ing time). We find that the global star formation 
history and the fraction of baryons converted into 
stars in each galaxy, /* = M*/Mb ar , is nearly inde- 
pendent on the assumed value of e* . This is typical 
for feedback regulated star formation. Hence, the 
star formation efficiency, /* , is not an assumed pa- 
rameter but it is extracted from the simulations. 
Thus, the derived star formation efficiency /* is a 
very generic prediction of our feedback model be- 
cause it is nearly independent of the assumed value 
of e*, that is instead quite uncertain. 

3. Contribution to Reionization. Due to the feedback- 
induced bursting mode of star formation in pre- 
reionization dwarfs, the cosmological H II regions 
that they produce remain confined in size and never 



reach the overlap phase {e.g., Fig. [T] above). Reion- 
ization is completed by more massive galaxies. 

4. Gas photo- evaporation and metallicity. Star- 
forming dwarf galaxies show large variations in 
their gas content because of the combined effects of 
stellar feedback from internal sources and IGM re- 
heating. Ratios of gas to dark matter lie below the 
cosmic mean in halos with masses Md m < 10 s M Q . 
Fig. |4jleft) shows the fraction of baryons retained 
by dark and luminous halos. It is clear that small 
mass luminous halos lose most of their gas before 
reionization due to internal radiation sources. Dark 
halos instead are able to retain gas for a longer 
time (see also iHoeft etldl I2006D . Gas depletion 
increases with decreasing redshift: the lower-mass 
halos lose all their gas first but, as the universe 
evolves, larger halos with Md m ~ 10 8 M Q also lose 
a large fraction of their gas. Gas photoevaporation 
reduces the metallicity of pre-reionization dwarfs to 
levels consistent with observations of dSph galax- 
ies. 

5. Number of luminous galaxies. Only about 1% of 
dwarf dark matter halos with mass M& m ~5x 10 6 
Mq, assembled prior to reionization, are able to 
form stars. The fraction of luminous halos scales 
with the halo mass as Fi um oc M| m . Thus, most 
halos with mass ^ 5 x 10 7 Mq are luminous. 
Fig. HJright) shows Fi um as a function of the halo 
mass at redshifts z — 12.5 and z = 10.2. The fig- 
ure also illustrated that /* depends on the environ- 
ment. Namely, it depends on the proximity of the 
pre-reionization dwarfs to other luminous galaxies. 
We find ~ 450 dwarf galaxies per Mpc 3 with bolo- 
metric luminosity between 10 4 and 10 s L . The 
luminosity function is rather flat at low luminosi- 
ties, with about 10 galaxies per Mpc 3 in the range 
10 7 < L < 10 8 L Q , and about 220 Mpc -3 in the 
range 10 4 < L < 10 5 L Q and 10 5 < L < 10 6 L . 
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Fig. 5. — (Left). Average star formation efficiency (i.e., fraction of the collapsed baryon mass converted into stars), (/*) = (M„ / Afbar) 
(thick curves), and gas fraction (i.e., fraction of the collapsed baryon mass retained in the gas phase) (/ g ) = (M g /M^ aI ) (thin curves), of 
the first galaxies as a function of their halo mass for run S2 in RGS08 at z = 14.5, 12.5, 10.2, and 8.3. For comparison, the symbol with error 
bar shows the expected star formation efficiency (roughly M» ~ 30 — 300 Mq divided by the baryonic mass of the halo) in the first mini 
halo of mass 10 6 M simulated by I Abel et all 1120031 . The dot-dashed line shows a power-law fit for the mean SFE, (/*)(M t j m ,z) oc M^ m . 
The SFE is nearly independent of redshift apart from an increase in halos with M^ m < 10 7 Mq at z ~ 8. (Right). Same as the left panel 
but showing the star formation efficiency, /*, for individual galaxies in the simulation at z = 10.2. Circles, from smaller to the larger, refer 
to galaxies with gas fractions / g < 0.1% (blue), 0.1% < / g < 1% (cyan), 1% < /g < 10% (red) and / g > 10% (green), respectively. 



6. Star formation efficiency and mass-to-light. The 
mean star formation efficiency (/*} = (M*/Mg^*), 
where M£^ x ~ Md m /7, is found to be nearly inde- 
pendent of redshift and to depend on total mass as 
(/*} oc Mj m with a = 2 if the radiative feedback 
is strong (i.e., top heavy IMF and/or large (fesc)) 
and a = 1.5 if the feedback is weak. This is shown 
in Fig. [Sfleft), where we plot the mean star forma- 
tion efficiency, (/*), and the mean gas fraction (f g ) 
in halos of mass M^ m . 

7. Scatter of the mass-to-light ratio. A tight rela- 
tionship between the star formation efficiency /* 
and the total mass of halos holds only for galax- 
ies with Afdm > 5 x 10 7 M Q . In lower-mass halos, 
the scatter around the mean (/*) is increasingly 
large (see Fig. Etright)). For a given halo mass, 
the galaxy can be without stars (dark galaxy) or 
have star formation efficiency /» ~ 0.1. However, 
only a few dark galaxies with mass at formation 
Afdm ~ 1 _ 5 x 10 7 M should exist in the Local 
Group. 

8. Size and morphology of stellar component. Galax- 
ies with masses Mdm < 10 8 Mg have a low sur- 
face brightness and extended stellar spheroid. At 
z ~ 10, the outer edges of the stellar spheroid 
nearly reaches the virial radius. In more massive 
galaxies that cool more efficiently by Lyman-alpha 
radiation, the stars and light are more centrally 
concentrated. Fig. [6] shows the structural proper- 
ties of the dark matter and stellar halo in three of 
the most massive galaxies in our simulation. These 
dwarf galaxies have properties similar to Draco and 
Umi dSphs. The figure also shows that the veloc- 
ity dispersion of the stars in these dwarfs is about 
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Fig. 6. — The top panels show the surface brightness (black point) 
and luminosity density (gray points) radial profiles for 3 galaxies 
from among the most massive dwarf galaxies extracted from the 
RG05 simulation at z = 8.3. We have evolved the stellar population 
passively to z = 0. The simulated galaxies shown in this figure 
have a stellar spheroid similar to Draco and Umi already in place 
at the time of formation (i.e., the spheroid is not produced by tidal 
effects). The bottom panels show the velocity dispersion profiles of 
the dark matter (gray points) and of the stars (black points) for the 
same 3 galaxies. The velocity dispersion of the stars is split in the 
radial and tangential components. All quantities are spherically 
averaged because the dark matter and stars have nearly spherical 
symmetry. 
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FlG. 7. — The evolution of the gas density in the core of 
dark halos due to cold accretion from the IGM for halos with 
D v ; r = 18,15,12,9 and 6 km s — 1 (from the top to the bottom), 
and corresponding to the circular velocities at the core radius 
^cir ('"core ) = 0.66f„i r (shown by the labels). The curves are assum- 
ing isothermal equation of state of the gas but the dashed portions 
show the parameter space in which such assumption fails because 
the gas cannot cool sufficiently fast as it is compressed toward the 
center of the halo. 

3.1.5. Photo- evaporation and Reionization feedback 

The small total mass of the first galaxies has two 
other implications. First, the ionizing radiation emit- 
ted by massive stars can blow out most of the gas be- 
fore SN-driven winds become important, further reduc- 
ing star formation rates (see RGS08). Second, the in- 
crease in temperature of the IGM to 10, 000 - 20, 000 K 
due to H I reionization prevents the gas from condensing 
into newly virialized halo s with circular velocities smaller 
than 10- 20 km s" 1 (e.o..lBabul & ReeslH9 92; Efst athiou l 



1992t iBullock et al 
200a IQkamoto et al.l 



20001: IGnedinl I2000bb iHoeft et all 



20081 ) . It follows that dwarf galax- 



< 10 — 20 km s lose most of their gas 
before reionization and stop accreting new gas and form- 
ing stars after reionization. 

The value i> max ~ 20 km s _1 that we use to define a 
fossil is motivated by the fundamental differences in cool- 
ing and feedback processes discussed above that regulate 
star formation in the early universe. It is not the critical 
value for suppression of gas accretion due to reionization. 
Indeed, we discuss in §|4]that pre-reionization fossils may 
have a late phase of gas accretion and star formation 
well after reionization, at redshift z < 1 — 2. Thus, a 
complete suppression of star formation after reionization 
(about 12 Gyr ago) is not the defining property of a fossil 
dwarf. 

4. LATE TIME COLD ACCRETION FROM THE IGM 

The ability of the IGM gas to condense at the center 
of dark halos depends on the ratio, T = w v ir/c s .igm, of the 
circular velocity to the IG M sound speed , and also on the 
dark halo concentration, c (jRicot ti 2009). Typically, the 
concentration of a halo is c v 



4 at the redshift o f viri- 
alization (jBullock et al.M200H : iWechsler et al.ll2002f) but, 
as the halo evolves in the expanding universe, its con- 
centration increases oc (1 + z v ir)/(l + z). The evolution 



of the halo concentration with redshift can be under- 
stood in the context of the theory of cosmolog ical sec- 
ondary infall of dark matter ( Bcrtschingcr 1985) and has 
been quantified using N-bo dy simulations (|Bullock et al.l 
!200lHWechsler et al.ll2002T) . Thus, primordial halos with 
u V ir < 10 — 20 km s" 1 stop accreting gas immediately 
after reionization, but, in virtue of their increasing con- 
centration and the decreasing temperature of the IGM at 
z < 3 (after He II reionization), they may start accr eting 
gas from the IGM at later times (see, iRicotti 2009]) . As 
a result, we expect that if the fossils of the first galax- 
ies exist in the Local Group (RG05), they may have a 
more complex star formation history than previously en- 
visioned. A signature of this model is a bimodal star 
formation history with an old (~ 13 Gyr) and a younger 
( ^ 5 — 10 Gyr, depending on the halo mass) population 
of stars. Leo T properties can be r eproduced by this 
simple model for late gas accretion (|Ricottil I2009D . In 
addition , Leo T seems to sho w a bimodal star formation 
history (jde Jong et al.l f2008h as expected in our model. 
Still, other models may als o explain the obser ved star 
formation history of Leo T ( Stinson et al.ll2007t ). 

Perhaps more important is the possible existence of 
dark galaxies: small mass halos containing only gas but 
no stars. Dark galaxies are most likely to exist if pre- 
reionization fossils do not form efficiently due to dom- 
inant negative feedback in the early universe (e.g., H2 
photo-dissociation by the FUV background). 

The late gas accretion from the IGM is shown in Fig- 
ure [7] for dark halos with circular velocity at virializa- 
tion w v ; r — 18,15,12,9 and 6 km s _1 (lines from the 
top to the bottom). The lines show the evolution of 
the gas density in the core of the halo as a function 
of redshift. The core radius is typically 100 pc and 
the labels show the circular velocity at the core radius 
(v C i T (r COIC ) « 0.66-Uot,. w 0.624u max , where v V ir is circular 
velocities at the virial radius and i> ma x is the maximum 
circular velocity) . We show the evolution of the gas 
density only for halos that are affected by reionization 
feedback. More massive halos will also accrete gas from 
the IGM as they evolve in isolation after virialization, 
but the gas accretion is continuous and not affected by 
reionization. 

Under the scenario in which halos with masses smaller 
than the critical value of 10 8 - 10 9 M remain dark due 
to feedback effects, the increase in their dark matter con- 
centration and the temperature evolution of the IGM will 
produce a late phase of gas accretion at redshift z < 1 — 2. 
If the gas has very low metallicity or is metal free, it is 
unlikely that the accreted gas will be able to form stars 
in the smallest mass halos. This is why we envisioned the 
possible existence of dark galaxies. However, their mass 
would be smaller than the mass of any luminous galaxy 
and their discovery would be challenging. 

The level of metal pre-enrichment necessary to ini- 
tiate star formation in minihalos that experience a 
late phase of gas accretion can be roughly estimated 
from the cooling function from hyperfine transitions 
of oxygen and carbon: A23 ~ 10~ 3 (Z/Zq), where 
A23 = 10~ 23 erg s" 1 cm 3 and Z is the gas metallic- 
ity. A necessary condition for star formation is t coo i « 
(0.7 yr) r/(n g . coro A_23) < ^h, that can be written as 
ng.corc > 0.03 cm- 3 (Z/10- 2 Z©) -1 . The left panel 
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Fig. 8. — (Left). The gas density, 

Tig, C orc (solid curves), and hydrogen column density, A^h — 27* c ?7.g iC orc (dashed curves), within the core, 
r c , of a minihalo at redshift z = as a function its circular velocity at r c . The minimum v c - lT in each curve is determined by the condition 
t rec /i H < 1, necessary for cooling to T gas ~ 10 4 K. The horizontal lines with arrows show the requirement for cooling to temperatures 
below 10 4 K, necessary for initiate star formation, for gas metallicity Z = 0.1 (lower line) and 0.01 Zq (higher line). (Right). Same as in 
the left panel but for minihalos whose gas is able to cool to T giCorc = 5000 K due to metal pre-enrichment (Z = 0.1 Zq). These halos are 
likely able to sustain a multi-phase ISM and form stars. 

in Figure [8] shows n gjCO rc and -/Vpt in minihalos that 
evolve isothermally at T ~ 10 4 K but that do not form 
stars (i.e., candidates for extragalactic CHVCs and dark 
galaxies). The horizontal lines show the requirement for 
metal cooling and star formation assuming gas metallic- 
ity Z = 0.1 and 0.01 Z©. The right panel in Figure H 



shows n g!Coro , TVjj and M^y n /M gaa (the dynamical mass 
to gas mass ratio) in the core of minihalos that are able 
to cool to Tg coio = 5000 K (roughly the temperature of 
the ISM in Leo T), and thus form stars. The symbols 
show the observed value for Leo T. 

5. COMPARISON OF THEORY AND OBSERVATIONS 

5.1. Number of fossils and non-fossil satellites in the 
Milky Way 

N-body simulations can be used to estimate the num- 
ber of dark halos in the Milky Way with maximum 
circular velocity w max > 20 km s _1 . However, there 
is a complication to this naive calculation. A signif- 
icant fraction of dark halos that today have v max < 
20 k m s" 1 were once mor e massive, due to tidal strip- 
ping (|Kravtsov et al.ll2004T ) . According to our definition, 
dwarf galaxies formed in these dark halos would not be 
pre-reionization fossils if they had at any time during 
their evolution v max (t) > 20 km s~ x (see § 13.1. 5| . If 
the number of observed Milky Way satellites exceeds the 
estimated number of these massive halos we must con- 
clude that at least a fraction of the observed Milky Way 
satellites are pre-reionization fossils. 

However, there is an assumption in this scenario. One 
must assume that the stars in these halos survive tidal 
stripping for as long as the dark matter. In this case 
tidally stripped halos may indeed account for a fraction 
or al l of the newly dis c overed ultra- faint dwarfs. How- 
ever, iPenarrubia et al.l (|2008l ) find that tidally stripped 



dark halos lose their stars more rapidly than they lose 
their dark matter. Thus, they may become dark ha- 
los even though they were initially luminous satellites. 
These dark halos should not be counted as ultra-faint 
dwarfs. 

Using res ults of published N-bod y simulations of the 
Milky Way, iBovill fc Ricottil (|2009l ) have estimated the 
number of dark halos that have or had in the past 
v m a.x(t) > 20 km s _1 (i.e., non pre-reionization fossils). 
In Table [T] we summarize the results of the counts for 
dark matter and luminous satellites for two large N-body 
simulations of a Milky Way typ e halo: the "Aquarius' 
simulation (iSprineel et al.l 12008) and the Via Lactea I 
simulation (Diem and et al.ll2007l ). 



TABLE 1 

Number of observed satellites versus number of dark 

HALOS WITH -Umax(t) > 20 KM S _1 (i.e., NON PRE-REIONIZATION 

fossils) for the Milky Way 



Distance 


Luminous 


Dark halos with v 


»i(t) > 20 km s" 1 


from 


dwarfs 


Via Lactea I sim. 


Aquarius sim. 


center 




today any time 


today any time 


< 200 kpc 


176 to 330 


14 36 ± 8 


34 91 ± 20 


< 417 kpc 


304 to 576 


28 73 ± 16 


69 182 ± 40 



The number of luminous satellites that exist within 
the Milky Way is highly uncertain bey ond a distance 
from the Galactic center of 200 kpc. iTollerud et alJ 
(2008), after applying incompleteness corrections, es- 
timated 304-576 satellites within 417 kpc and about 
176-330 within 200 kpc (the numbers are from their Ta- 
ble 3). As shown in Table [TJ the existence of some pre- 
reionization fossils among the ultra-faint dwarfs appears 
to be favored by the data. However, the current uncer- 
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tainties on the completeness corrections of observations 
and on the simulations are too large to deem the exis- 
tence of fossils as necessary. 

The error bars on the theoretical estimate of the num- 
ber of fossils in the Milky Way shown in Table [1] come 
from uncertainties in the fraction of halos that were more 
massive in t he past. This fraction was derived from sim- 
ulations by (|Kravtsov et al.ir2004f ). Another uncertainty 
in the simulation results can be attributed to the differ- 
ent predictions for the number of Milky Way satellites in 
the Via Lactea I and II and Aquarius simulations. The 
discrepancy can be partially attributed to different cos- 
mology in the simulations but mostly because the Via 
Lactea I simul ation likely used erro neous initial condi- 
tions. Finally, I Toller ud et al.l ([2008) corrections on the 
number of observed satellites also rely on the radial dis- 
tribution of dark matter sub halos extracted from Via 
Lactea I simulations that may be erroneous. Once the 
discrepancies among different simulations are better un- 
derstood the number of simulated satellites of the Milky 
Way may be known with greater certainty. 

Using comparisons between the predicted and observed 
Galactocentric di stributions of dwarf satell ites around 
the Milky Way, iGnedin fc Kravtsovl (|2006l ). hereafter 
GK06, have estimated that pre-reionization fossils may 
constitute about 1/3 of Milky Way dwarfs. GK06 esti- 
mated the number of fossils in the Milky Way using data 
from the simulations of the first galaxies in RG05. GK06 
defined a fossil as a simulated halo which survives at 
z = and remains below the critical circular velocity of 
20 km s _1 with no appreciable tidal stripping (the usual 
definition of fossil adopted in this paper as well). They 
calculate the probability, P$ (w max , r), of a luminous halo 
with a given maximum circular velocity w max to survive 
from z = 8 (the final redshift of the RG05 simulation) 
to z — 0. The surviving halos are assigned a luminosity 
based on the Ly versus v max relationship from RG05. At 
z = 0, GK06 has a population of dwarf galaxies with a 
resolution limit of u max = 13 km s _1 . This limit cor- 
responds to a lower luminosity limit of Ly ~ 10 5 Lq, 
which includes Leo T and Canes Venatici I, but excludes 
all the other new ultra-faint Milky Way satellites. 

In Figure [5J we show the cumulative luminosity func- 
tion from GK06 for the Milky Way and M31 satellites. 
The lower panel shows satellites with distance from their 
host d < 100 kpc, the middle panel d < 300 kpc and the 
upper panel d < 1 Mpc. The gray lines show the GK06 
predictions, and the shaded region encompasses the er- 
ror bars. The resolution limits in GK06 cause halos with 
w mal < 17 km s _1 to be preferentially destroyed by tidal 
effects. The dashed line show the predicted luminosity 
function corrected for the resolution effects. Both the 
uncorrected (solid lines) and corrected (dashed lines) lu- 
minosity functions are plotted in the lower panel. The 
points with error bars show the observed luminosity func- 
tion of dSph galaxies around the Milky Way and M31 
corrected only for limits in sky coverage of the SDSS 
survey. The plot is from GK06 but has been updated to 
include the new ultra- faint dwarfs with Ly ^ 10 L©. 

The results of this model are consistent with the obser- 
vations. The model reproduces the Galactocentric distri- 
bution of the most luminous dSphs, even though in this 
model dSphs are not tidally stripped dlrrs. It also shows 
a good agreement with observations for luminosities that 
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Fig. 9. — Luminosity function of pre-reionization fossil dwarfs 
predicted in GK06 (red bands) plotted with the luminosity function 
for Local Group dSphs (points with error bars). The data from 
observations are corrected only for limits in sky coverage of the 
SDSS survey. 

can be considered nearly complete within a given Galac- 
tocentric distance. 

5.2. Statistical Properties of pre-reionization Fossils 

In this section, we compare the properties of the new 
dwarf galaxies discovered in the Local Group to the the- 
oretical predictions of simulations of primordial galaxies 
formed before reionization. The argument that justifies 
this comparison is that star formation stops or is greatly 
reduced after reionization (but see §2]). We do not expect 
two perfectly distinct populations of fossil galaxies with 
Umax < 20 km s" 1 and non-fossils with v max > 20 km s _1 , 
but a gradual transition of properties from one popula- 
tion to the other. Some fossils may become more massive 
than i> max ~ 20 km s _1 after reionization, accrete some 
gas from the IGM, and form a younger stellar popula- 
tion. These dwarfs are no longer defined as "fossils" . 
However, if the dark halo circular velocity remains close 
to 20 km s _1 the young stellar population is likely to be 
small with respect to the old one. In RG05 we call these 
galaxies "polluted fossils" because they have the same 
basic properties of "fossils" with a sub-dominant young 
stellar population. A similar argument can be made re- 
garding the late phase of gas accretion that may produce 
objects similar to Leo T. 

In Figs. nUirm we compare the RG05 predictions for 
the fossils of primordial galaxies to the observed proper- 
ties of the new Milky Way and M31 dwarfs. The sym- 
bols and lines in Figs. ITDlfTTI have the following meanings. 
All known Milky Way dSphs are shown by circles; An- 
dromeda's dSphs satellites are shown by triangles; sim- 
ulated fossils are shown by the small solid squares. The 
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Fig. 10. — (Left). Surface brightness and core radius vs. V-band luminosities. Small filled squares are simulated pre-reionization fossils 
from RG05, asterisks are dlrrs, crosses are dEs, closed circles are the previously known dSph around the Milky Way, closed triangles are 
previously known dSph around M31, and open circles and triangles are new dSph around the Milky Way an d M31 respectively. T he solid 
lines roughly show the detection limits inherent to the methods used to find the ultra-faints in the SDSS data (Koposov ct al. 2008) and the 
dashe d lines show the scaling relationships for more luminous Sc-Im galaxies (10 8 Lq ,$ Lb ~ 10 1:L Lp)) derived by [Kormcndy & F reeman! 
(2004). (Right). Mass-to-light ratio and velocity dispersion of a subset of the new dwarfs (Martin ct al. 2007; Simon L denall2007l) vs 
V-band luminosity. The symbols and lines are as in the left panel. 



solid and open symbols refer to previously known and 
new dSphs, respectively. The transition between fossils 
and non-fossil galaxies is gradual. In order to illustrate 
the different statistical trends of "non-fossil" galaxies we 
show dwarf irregulars (dlrr) as asterisks and the dwarf 
ellipticals (dEs) as crosses, and we show the statistical 
trends for more luminous galaxies as thick dashed lines 
on the right side of each panel. 

Figure [TU^left) shows how the surface brightness (top 
panel) and half light radius (bottom panel) of all known 
Milky Way and Andromeda satellites as a function of V- 
band luminosity compares to the simulated fossils. The 
surface brightness limit of the SDSS is shown by the thin 
solid lines in both panels of the figure. The new dwarfs 
agree with the predictions up to this threshold, suggest- 
ing the possible existence of an undetected population of 
dwarfs with Sy below the SDSS sensitivity limit. The 
new M31 satellites have properties similar to their previ- 
ously known Milky Way counterparts (e.g., Ursa Minor 
and Draco). Given the similar host masses and environ- 
ments, further assuming similar formation histories for 
the halos of M31 and the Milky Way, we may be tempted 
to speculate on the existence of an undiscovered popula- 
tion of dwarfs orbiting M31 equivalent to the new SDSS 
dwarfs. 

The large mass outflows due to photo-heating by mas- 
sive stars and the subsequent suppression of star forma- 
tion after an initial burst make reionization fossils among 
the most dark matter dominated objects in the universe, 
with predicted mass-to-light ratios as high as 10 4 and 
Ly ~ 10 3 — 10 4 L Q . Figure [TOMght) shows the veloc- 
ity dispersion (bottom panel) and mass-to-light ratios, 
M CT /Lv (top panel), as a function of V-band luminosity 
of the new and old dwarfs from observations in compari- 
son to simulated fossils. The symbols are the same as in 
the previous figures. Theoretical and observed dynam- 



ical masses are calculated from the velocity dispersions 
of stars (i.e., M a — 2r 1 / 2 & 2 /G), and do not necessarily 
reflect the total mass of the dark halo at virialization. 

Observations show that the value of the dynamical 
mass within the stellar spheroid, M ~ (1 ± 5) X 1 7 M^, 
rema ins relatively cons tant as a function o f Ly ([Mated 
119981) . Recent work bv lStrigari et all (120081) show s anal- 
ogous results to the one found bv lMateor(ri99l ). The 
dynamical mass of dwarf spheroidals within a radius of 
300 pc is relatively constant: M ~ 10 7 M©. The radii 
of the stellar spheroids in these dwarf galaxies may be 
either larger or smaller than 300 pc. In the later case, 
the determination of the mass of the dwarfs is uncertain. 

Our simulation provides some insight into the reason 
why the dynamical mass remains roughly constant in 
dSphs. The simulations show that in pre-reionization 
dwarfs, the ratio of the radius of the stellar spheroid 
to the virial radius of the dark halo decreases with in- 
creasing dark halo mass (i.e., the stellar profile becomes 
more concentrated for more luminous dwarfs). Thus, 
as the halo mass and virial radius increases, the stel- 
lar spheroid becomes increasingly concentrated in the 
deepest part of the potential well. If follows that the 
ratio, f a = M ff /M(i m , of the dynamical mass within the 
largest stellar orbits to total dark matter mass is also re- 
duced. Thus, the decrease of f a for increasing dark mat- 
ter mass of halos maintains the value of the dynamical 
mass M a = f a Mdm (measured by the velocity dispersion 
of the stars) almost constant, even though the total mass 
of the halo increases. The extent of the stellar spheroids 
in the lowest mass dwarfs is comparable in size to their 
virial radii at formation (see § 13. ip . 

The metallicity-luminosity relation of the observed and 
simulated dwarfs is shown in the left panel of Figure 111! 
[Fe/H] is plotted against V-band luminosity in solar 
units. Symbols for the previously known dwarfs, the new, 
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ultra-faint dwarfs, and simulated fossils are the same as 
in Figure [TUJ In this plot we color code simulated fos- 
sils according to their star formation efficiency, /*. Red 
symbols show simulated dwarfs with /» < 0.003, blue 
0.003 < /* < 0.03 and green /* > 0.03. 

Using the d ata for the metallicity collected in 

iBovill fc RicQttH |2009), the new ultra-faint dwarfs do 
not appear to follow the tight luminosity-metallicity re- 
lationship observed in more luminous galaxies (although 
error bars are large). Note that here, as well as in 
IBovill fc Ricottil ((2009) (although Table 3 in that paper 
was erroneously no t updated), we have plotted data from 
iKirbv et al] (|2008f) for the 6 ultra-faint Milky Way satel- 
lites provided in that paper. There are several physi- 
cal mechanisms that may produce the observed scatter 
in metallicities of dwarfs at a given constant luminosity. 
The large spread of star formation efficiencies produc- 
ing a dwarf of a given luminosity in our simulations is 
responsible for at least part of the large spread of the 
luminosity-metallicity relation. At this point it is un- 
clear whether our simulations can reproduce the scatter 
of metallicities of simulated fossils, or if perhaps the lu- 
minosity of the lowest luminosity ultra-faint dwarfs has 
been reduced due to tidal interactions. As mentioned be- 
fore we have suggested that the lowest luminosity ultra- 
faint dwarfs have not yet been discovered because their 
surface brightness lies below the SDSS detection limits. 

Figure [TIT right) shows the scatter of the metallicity of 
the stars, crpe/H] > plotted against V-band luminosity and 
[Fe/H] respectively. The various point types and colors 
are the same used in the left panel of figure[llj The large 
spread in the metallicity of the stars is a natural conse- 
quence of the hierarchical assembly of dwarf galaxies in 
cosmological simulations. It is not necessarily an indica- 
tion that star formation in dwarf sat ellites was protracted 
for lon ger than 1 Gyr, as argued in lGrebel fc Gallagherl 
(2004) to prove that star formation in dwarf spheroidals 
is not stopped by reionization feedback. 

6. DISCUSSION: THE TIDAL SCENARIO VS FOSSIL 
HYPOTHESIS 

According to the results summarized above in Table [T] 
the number of dark matter satellites of the Milky Way 
that have or had in the past v max > 20 km s _1 is smaller 
than the number of observed luminous satellites (after 
applying completeness corrections). This implies that 
non-fossil galaxies cannot account for all the observed 
satellites. Thus, pre-reionization fossils should exist. 

However, we have already discussed the several uncer- 
tainties in estimating the numbers summarized in Ta- 
ble [1] Additional uncertainties that render the identi- 
fication of fossils uncertain are the following. The mass 
and circular velocity of the dark matter halo of the Milky 
Way are not known precisely. Simulations should take 
into account the effect of baryons in modifying the den- 
sity profile and the radial distribution of satellites. The 
effect of tidal stripping on the properties of the stars 
in the satellites is not well understood, thus we do not 
know if the tidal scenario is consistent with observations 
of ultra- faint dwarfs. The luminosity and stellar proper- 
ties of non-fossil dwarf satellites are not known. 

Non- fossil galaxies with v max > 20 km s _1 may lose 
a substantial fraction of their mass due to tidal inter- 
actions. If they survive the interaction, their proper- 



ties, such as surface brightness and half light radius, may 
be modified. iKravtsov et al.1 (|2004f) estimate that about 
10% of Milky Way dark matter satellites were at least 
ten times more massive at their formation than they are 
today. Although their simulation does not include stars, 
they favor the idea that the stellar properties of these ha- 
los would remain unchanged (i.e., dSphs are not tidally 
stripped dlrr). In their model the majority of brighter 
dwarf satellites have been considerably more massive in 
the past and could have formed their stars undisturbed 
after reionization. More precisely, the redshift of reion- 
ization does not affect the results of their model for clas- 
sic dwarfs because the probability of these to be fossils 
is low. 

This version of the tidal model may be hard to dis- 
tinguish observationally from the model we propose for 
the fossils because in both models the properties of the 
dwarfs are not modified by tidal forces (i.e., their prop- 
erties are those at formation). In addition, fossil dwarfs 
may stop forming stars for only about 2 to 4 Gyr after 
reionization, before starting to accrete gas again from the 
IGM. Thus, reionization may imprint a bimodal star for- 
mation history in some fossils, but this signature is not 
a robust discriminant because the star formation history 
of dwarfs cannot be determined with sufficient accuracy. 

Observations seem to suggest that star formation in 
dwarf galaxies slightly more massive than 10 8 — 10 9 M Q 
may be similar to star formation in fossils and thus fit 
the observed properties of classical dSphs without re- 
quiring significant ti dal stripping of s tars. If star forma- 
tion was included in IKravtsov et alj (|2004h . their model 
may have reproduced the properties of observed dwarf 
satellites that our simulations of pre-reionization dwarf 
galaxies already does. The differences between the two 
models will depend on whether fossil galaxies are allowed 
to form and on their properties. In other words, the two 
models may differ on the assumed mass of the smallest 
dark halo that can host luminous satellites. This critical 
mass cannot be directly observed in dwarf galaxies but, in 
principle, can be constrained by comparing the observed 
number of luminous satellites to the model predictions. 
Determining the minimum mass for a dark halo to be- 
come luminous is of great importance in understanding 
galaxy formation in the early universe. 

To summarize, there are a few observational tests that 
can be used to distinguish true fossils from dSphs or dEs 
that form in more massive halos and form stars unaf- 
fected by IGM reionization. True fossils should have ei- 
ther a single old stellar population or have a bimodal star 
formation history produced by a temporary suppression 
of star formation after reionization and late gas accre- 
tion. In addition, if the number of observed Milky Way 
satellites (or the number of isolated dwarfs) exceeds some 
critical value determined using N-body simulations (e.g., 
see Table [T]), we may conclude that some pre-reionization 
fossils do exist in the Local Group. 

It is likely that these tests will prove inconclusive for 
some time to come, unless the number (after corrections 
for completeness) of new ultra-faint galaxies surges in the 
coming years. The weakness of the star formation his- 
tory test is that it requires measurements with precision 
of 1-2 Gyr of the stellar populations in order to be really 
discriminating between models that are quite similar to 
each other. This is hard to achieve especially for ultra- 
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Fig. 11. — (Left). Mean metallicity of the stars vs. V-band luminosity for Local Group dSphs plotted against RG05 predictions. They 
symbols have the same meaning as in Fig. 8. In this plot simulated fossils (shown as small squares) are color coded according to their star 
formation efficiency, /*: red symbols show simulated dwarfs with /„ < 0.003, blue 0.003 < /* < 0.03 and green /„ > 0.03. (Right). Same 
as the left panel but showing the spread of the stellar metallicities in each dwarf (i.e., variance of the metallicity distribution) vs. their 
V-band luminosity. 



faint dwarfs with few stars. If the number of ultra-faint 
dwarfs remains about the same as today, the number ar- 
gument may also remain controversial until more detailed 
theoretical modeling can reduce the current uncertain- 
ties surrounding the expected number of dark halos in 
the Milky Way and the completeness corrections of the 
observations. Ultimately, the case for the origin of ultra- 
faint dwarfs must be made on the basis of the model that 
does the best job of reproducing available observations. 

Finally, even if pre-reionization fossils do not exist 
(i.e., halos with v max < 20 km s _1 are all dark), a frac- 
tion of them should be able to accrete some gas at red- 
shift z < 1 — 2 and might be discoverable in the outer 
parts of the Local Gro up using Ha or 21 cm sur veys 
{e.g., ALFALFA survey. lOiovanelli et aT1[2fM 120071 ) . Of 
course, one should prove that the gas clouds are embed- 
ded in dark halos. Measurements of the gas cloud size, 
column density, and velocity broadening of the emis- 
sion/absorption lines can be used to discriminate be- 
tween "dark galaxies" and tidal debris. This is because 
the gas in dark galaxies is confined by the gravitational 
potential of the dark matter halo, while tidal debris or 
clouds formed via ther mal instability are confined by the 
external gas pressure ()Ricottil [20091. This is another 
promising direction for determining the minimum mass 
of luminous galaxies in the universe. 

Another variant of the tidal hypothesis for the origin of 
dSphs is a scenario in which dlrr galaxies transform into 
dSphs as they fall into the Milky Way and Andromeda, 
due to tid al and ram pressure stripping dMaver et al.l 
120071120061 ). A work bv lPeiiarrubia etaTI (|2008D explores 
the idea that ultra-faint dSphs are tidally stripped dlrrs. 
They achieve some success in reproducing observed prop- 
erties of ultra-faint dwarfs. While this type of tidal strip- 
ping can reproduce properties of an individual galaxy, it 
is unable to completely reproduce all the trends in the 
ultra-faint population. This is primarily seen in the kine- 
matics of the ultra- faint dwarfs. Tidal stripping predicts 



a steeper than observed drop i n the velocity dispersio n 
of the stars with decreasing Ly (jPenarrubia et aLl l2008) . 
In addition several dSph do not show signs of strong 
tidal stripping. And XII and And X IV may be on their 
first approach to the Local Group (|Martin et al.l [20061 : 
iChapman et al.l I2007T ) ■ Other examples of dSphs that 
are found distant from the center o f their host galaxies 
are A nd XVIII, Cetus and Tucana (jMcConnachie et alJ 

Hoop). 

Finally, another interesting case study is Leo T, that 
we have discussed extensively above in § 12.21 and § 2J 
Leo T properties can be explained in some detail as being 
a fossil that experienced a late phase of gas accretion 
(|Ricottill2009| ). However, another possibility that should 
be explored quantitatively with simulations is that Leo T 
is more massive than a fossil but less massive than dlrr 
galaxies. 

7. CONCLUSIONS AND FUTURE WORK 

We have summarized our work on the formation of the 
first galaxies before reionization (i.e., pre-reionization 
dwarfs) and the quest to identify the fossils of these 
first galaxies in the Local Group. The definition of a 
pre-reionization fossil is not directly related to the sup- 
pression of star formation experienced by these galax- 
ies due to reionization feedback. Indeed, we discussed 
how pre-reionization fossils may experience a late phase 
of gas accretion and possibly star formation at redshift 
z < 1 — 2. Most importantly, fossils are a population of 
dwarf galaxies whose formation (i.e., the fraction of ha- 
los that are luminous) is self-regulated on cosmological 
distance scales by radiative processes. Their existence 
is not certain due to a possible strong negative feedback 
that may prevent the majority of these halos from ever 
forming stars. In addition, if negative feedback heav- 
ily suppresses the number and luminosity of these first 
galaxies, more massive halos with u max > 20 km s^ 1 will 
evolve differently because of the lower level of metal pre- 
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enrichment of the IGM. To summarize, the critical circu- 
lar velocity u max ~ 20 km s -1 that we adopt to define a 
fossil is primarily motivated by fundamental differences 
in cooling and feedback processes that regulate star for- 
mation in these halos in the early universe. However, it 
is also close to the critical va lue for continu e d gas accre- 
tion after IGM reionization (|Gnedinll2000bl : iHoeft et all 
l200l lOkamoto et al.ll2008h . 

The number of Milky Way and M31 satellites provides 
an indirect test of galaxy formation and the importance 
of positive and negative feedback in the early universe. 
This test, although the uncertainties are large, supports 
the idea that a fraction of the new ultra-faint dwarfs 
are fossils. The good agreement of the SDSS and new 
M31 ultra-faint dwarf properties with predictions of our 
simulations (RG05, GK06, Bovill & Ricotti 2009) does 
not prove the primordial origin of the new ultra-faint 
dwarfs, but it supports this possibility. 

More theoretical work and more observational data are 
needed to prove that some dwarfs in the Local Group are 
true fossils of the first galaxies. Future theoretical work 
should focus on improving the accuracy of predictions on 
the properties of dwarf galaxies formed before reioniza- 
tion and their evolution to the present day. Modeling the 
evolution of the baryonic component after reionization in 
dwarf satellites and in the Milky Way - Andromeda sys- 
tem may be necessary to make robust predictions. More 
observational data will certainly be available in the near 
future. A large number of surveys, both at optical and 
radio wavelengths will be online in the near future {e.g., 
Pan-STARRS, LSST, ALMA, EVLA, JWST, SKA to 
mention a few). Different survey strategies may be used 
to find and characterize fossil dwarf galaxies. A deep pen- 
cil beam survey would be useful to find the faintest dwarf 
satellites of the Milky Way and determine more precisely 
their Galactocentric distribution. A shallower all sky sur- 



vey could be used to quantify the degree of anisotropy in 
the distribution of satellites around the Milky Way. 

The star formation history of the dwarf galaxies is not 
strongly discriminatory because fossil galaxies may have 
a late phase of gas a ccretion and star formation during 
the last 9-10 Gyrs (|Ricottill2009h . The distinction be- 
tween fossils and non-fossil galaxies may be quite elusive 
but it is nevertheless important to understand galaxy for- 
mation and feedback in the early universe. Arguments 
based on counting the number of dwarfs in the Local uni- 
verse are among the more solid arguments that could be 
used to prove the existence of fossil galaxies (see Table [T]). 

Future tests may be provided by deep surveys looking 
for ultra-faint galaxies in the local voids or looking for 
gas in dark galaxies {i.e., dark halos that have been able 
to accrete gas from the IGM at z < 1 — 2). Ultra- faint 
dwarfs should be present in the voids if dwarf galaxies 
formed in large numbers before reionization (Bovill & 
Ricotti, in preparation). If pre-reionization dwarfs never 
formed due to dominant negative feedback in the early 
universe, it is possible that a faint (in Ha and 21 cm 
emission) population of dark galaxies exists in the outer 
parts of the Local Group. Hence, another way to detect 
fossil galaxies in the outer parts of the Milky Way or 
outside the super-galactic plane would be to search for 
neutral or ionized gas that they may have accreted from 
the IGM. Future radio telescopes {e.g., ALMA, EVLA, 
SKA) may be able to detect neutral hydrogen in dark 
galaxies or in ultra- faint dwarfs. Ionized gas in the outer 
parts of dark halos may be observed in absorption along 
the line of sight of distant quasars (for instance in O VI 
or O IV with COS on the HST). However, the probability 
that a line of sight toward a quasar intersects the ionized 
gas collected from the IGM by dark or fossil galaxies 
might be small. Additional theoretical work is required 
to address these issues. 
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